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Abstract 



• We explore the use of polarized e^/e beams and/or the information on final state 

decay lepton polarizations in probing the interaction of the Higgs boson with a pair of 
Oh| vector bosons. A model independent analysis of the process e^e~ — > ffH, where / is 

any light fermion, is carried out through the construction of observables having identical 
properties under the discrete symmetry transformations as different individual anomalous 
^ interactions. This allows us to probe an individual anomalous term independent of the 

^ ' others. We find that initial state beam polarization can significantly improve the sen- 

. sitivity to CP-odd couplings of the Z boson with the Higgs boson (ZZH). Moreover, 

I an ability to isolate events with a particular r helicity, with even 40% efficiency, can im- 

' prove sensitivities to certain ZZH couplings by as much as a factor of 3. In addition, 

. the contamination from the ZZH vertex contributions present in the measurement of 

O I the trilinear Higgs-W^ (WWH) couplings can be reduced to a great extent by employing 

polarised beams. The effects of initial state radiation (ISR) and beamstrahlung, which 
can be relevant for higher values of the beam energy are also included in the analysis. 



^ . 1 Introduction 

Although the Standard Model (SM) has been highly successful in describing all the available 
experimental data, the precise mechanism of the breaking of the SU(2) ®U(1) gauge symmetry 
and consequent generation of masses for all the elementary particles is still very much an open 
question. In the SM, the symmetry is broken spontaneously giving rise to masses for all the 
elementary particles via the Higgs mechanism thereby requiring the presence of a spin-0 CP- 
even particle, namely the Higgs boson [1-4] . However, so far, there exists no direct experimental 
evidence for the the same. Not surprisingly, therefore, the search for the Higgs boson and the 
study of its various properties, comprise one of the major aims of all the current and future 
colliders [5]. The Large Hadron Collider (LHC), soon to go into operation, is expected to shed 
light on the mechanism of electroweak symmetry breaking (EWSB). It is designed to be capable 
of finding the SM Higgs boson over most of the theoretically allowed range for its mass [6]. 

Direct searches at the LEP gives a lower bound on the mass of the SM Higgs boson: run > 
114.4 GeV [7]. On the other hand, electroweak precision measurements put an upper bound 
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on its mass of about 182 GeV at the 95% confidence level (CL) [8]. Note, though, that these 
mass bounds are model-dependent and various extensions of the SM admit different allowed 
ranges of mn- For example, the lower bound can be relaxed in generic two-Higgs doublet 
models (2HDM) [9] and more spectacularly in multi-Higgs models with CP violation [10,11]. 
As a matter of fact, even in the Minimal Supcrsymmetric extension of the Standard Model 
(MSSM) [12], once additional CP- violation is admitted in the scalar sector [10], direct searches 
at LEP and elsewhere still allow a Higgs boson mass as low as 10 GcV [13]. Similarly, the non- 
minimal supcrsymmetric standard model admits very light spin-0 states even without invoking 
additional sources of CP violation [14]. In certain extensions, the upper bound on the mass of 
the (hghtest) Higgs boson may also be substantially higher [15]. A more detailed discussion of 
the subject may be found in Refs. [4, 14]. 

While the SM contains only a single CP-even scalar state, in general, various extensions of 
the SM mentioned above contain more than one Higgs boson and some, possibly, with different 
CP properties. For example, the 2HDM — of which the MSSM is a particular case — consists 
of five spin zero particles: two CP-even neutrals, a CP-odd neutral and a pair of charged 
scalars. If the MSSM parameters admit CP violation, the neutral particles may no longer be 
CP eigenstates. The aforementioned dilution of the experimental lower limits is, generically, 
the result of a reduced coupling of the lightest spin-0 state with the Z due to the mixing of the 
SM higgs with the other (pseudo) scalars in the model. 

Thus, even after the LHC sees a signal for a Higgs boson, a study of its properties (including 
CP) and precise measurements of its interactions would be necessary to establish the nature 
of electrowcak symmetry breaking. Such a detailed study of this sector may also provide 
the footprints of new physics beyond the SM. This, though, would be possible only at the 
International e~^e~ Linear Collider (ILC) [16, 17] and combined information from the ILC and 
LHC [18, 19], will be required to estabhsh it as the SM Higgs boson. A key step in this direction 
is the determination of the tensor structure of the coupling of the spin-0 state with the different 
SM particles. A model independent analysis would, then, incorporate the most general form 
for this tensor structure as allowed by symmetry principles, the anomalous parts having been 
assumed to have come from effects of high scale physics. The couplings of the Higgs boson with 
a pair of gauge bosons V{V = 7, 1^ and Z) as well as that with a tt pair have been studied very 
thoroughly in the context. The tensor structure can be inferred from kinematical distributions 
and polarisation measurements for various final state particles. In this study, we concentrate 
on the trilinear Higgs-Z {ZZH) and Higgs- (WWH) coupling, in particular focusing on the 
utility of beam polarisation, measurement of the final state particle polarisation as well as the 
use of higher beam energies, for the process e+e" ffH. 

At an e'^e" collider, the Z boson produced in the process e'^e" — > ZH is, at high energies, 
longitudinally polarised when produced in association with a CP-even Higgs boson and trans- 
versely polarised in case of a CP-odd Higgs boson. The angular and energy distribution of the Z 
boson can, thus, provide a wealth of information about the ZZH coupling [20,21]. The shapes 
of the threshold excitation curve in the processes e+e~ — >■ ZH [22,23] and e^e~ — > ttH [24] 
constitute model independent probes of the tensor structure of the ZZH and the ttH cou- 
pling respectively. Many detailed studies of how kinematical distributions for the processes 
e~^e~ — > ffH, proceeding via vector boson fusion and Higgs-strahlung can be used to probe 
the ZZH vertex exist [25,26]. The anomalous ZZH vertex, in the context of higher dimen- 
sional operators has been studied in Refs. [26-35] for a Linear Collider (LC). Ref. [31] is one 
of the pioneering studies and contains a very extensive analysis, using the optimal observable 
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technique [36], to probe ZZH and ^ZH couplings, whereas Refs. [32-34] use asymmetries 
constructed using differences in the kinematical distributions of the decay products. 

The ZZH vertex could be probed at the LHC in a similar fashion, again using kinematic 
distributions, threshold behaviour as well as asymmetries in the Higgs decays [37-43]. This, 
alongwith the WWH vertex, can also be studied through vector boson fusion at the LHC [44,45]. 
Angular distributions of the decay products have been used in Ref. [46] to study the VVH 
(where V = Z/W) vertex in the process j'-f —>■ H —>■ W^W^ / ZZ. 

In Ref. [34], an exhaustive set of asymmetries, which could probe each of the ZZH anoma- 
lous coupling independent of the others, were constructed. Defining kinematical observables 
which are either odd or even under the different discrete symmetry transformations, the said 
asymmetries are the expectation values of the sign of these observables. In the approxima- 
tion of small contribution from anomalous parts (which amounts to retaining terms only upto 
the linear order in the anomalous couplings), these asymmetries are then proportional to the 
coefficient of the term in the Lagrangian with the corresponding transformation properties. 
However, many of these asymmetries turned out to be proportional to the difference between 
the squared right and left handed couplings of the fermion to the Z boson (viz. — r^), and 
consequently were rather small, on account of the electrons (charged leptons) being involved. 
It follows then that the sensitivity of these asymmetries to the anomalous couplings could be 
enhanced by either using the polarized beams or through a measurement of the polarization of 
the final state particles. 

In the unpolarized case, the determination of the anomalous WWH coupling suffers a 
large contamination from the contribution from the s-channel diagram, arising from the ZZH 
coupling. We look at the possibility of reducing this contamination by the use of polarised 
beams. For completely polarized and e~ beams, (Tlr receives contributions from both 
the Bjorken (s-channel) and fusion (t-channel) diagrams, whereas only the s-channel diagram 
contributes to a^L. Thus, beam polarisation may also be used to enhance the sensitivity to 
WWH anomalous couplings, and this constitutes part of our investigations. Furthermore, we 
also study the dependence of the sensitivities on the beam energy; again with an aim to enhance 
the t channel contribution and hence the sensitivity to the WWH couplings. Effects of both 
initial state radiation (ISR) [47] and beamstrahlung [48,49] have been included in this study as 
they ought to be. 

The rest of the paper is organized as follows: in Sec. [2] we discuss possible sources of 
anomalous VVH couplings. Various kinematical cuts on final state particles used to suppress 
the background are discussed in Sec. [31 The ZZH vertex is examined in detail in Sec. HJ 
with the various observables being defined in Sec. 14. the effects of beam polarization being 
discussed in Sec. 14.21 and the use of final state r-polarization in Sec. 14. 3[ In Sec. 14.41 we discuss 
the improvements possible in the reach for the anomalous ZZH coupling using final state r 
measurement with polarised initial beams. In Sec. [5] we construct some observables using the 
polarisation of initial beams to constrain the WWH couplings. In Sec. |6] we present results 
on the dependence of the sensitivities to different couplings to the beam energy, including the 
effects of ISR and beamstrahlung. Finally we summarize our results in Sec. [71 
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2 The VVH Couplings 



Within the SM/MSSM, the only interaction term involving the Higgs boson and a pair of gauge 
bosons arises from the Higgs kinetic term in the Lagrangian. However, once we accept the SM 
to be only an effective low-energy theory, higher-dimensional (and hence non-renormalizable) 
terms are allowed. The most general VVH vertex, consistent with Lorentz invariance and 
current conservatioiJl] can be written as 



^uu — gv 



O'V Qiiu H ^-n-{kiyk2a — Qav ^1 " ^2) H V ^l^yaP ^"^2 

my my 



where ki denote the momenta of the two VT's (2''s). Here 

= e cot OwMz , gf" = 2 eMz / sin 2ew , 

9w being the weak- mixing angle and ef^ua/s the antisymmetric tensor with 60123 = 1- Within 
the SM, az = aw = 1 and by = by = at tree level. Anomalous parts may arise on account of 
higher order contributions in a renormalizable theory [50] or from higher dimensional operators 
in an effective theory [51]. While the imposition of SU(2)l U(1)y invariance relates the WWH 
couplings to the ZZH ones, such an assumption restricts the nature of the physics beyond the 
SM. Instead, we view them purely as phenomenological inputs and study their effect on various 
final state observables in collider processes. 

In general, each of these couplings can be complex, reflecting final state interactions, or 
equivalently, absorptive parts of the loops either within the SM or from some new physics beyond 
the SM. However, for each of the observables that we construct for the process e~^e~ — >■ ffH, 
one overall phase can always be rotated away and we may choose that to be corresponding to 
either or aw In our analysis, we choose az to be real and allow the others to be complex. 
Further we also assume az and a^ to be close to there SM value i.e. ay = 1 + Aay- 

For a generic multi-doublet model, supersymmetric or otherwise, couplings of the neutral 
Higgs bosons to a pair of vector bosons (V = Z, W) obey a sum rule [52-54]: 

E4yH, = 1- (2) 

i 

Although avvHi for a given Higgs boson in different models, such as MSSM, can be significantly 
smaller than the SM value, the presence of higher SU{2)l multiplets or more complicated 
symmetry breaking structures (such as those within higher-dimensional theories) [54] would 
lead to more complicated sum rules. 

The terms containing ay and by in Eq. [1] constitute the most general coupling of a CP- 
even Higgs boson with two vector bosons whereas the 6y term corresponds to the CP-odd one. 
Simultaneous presence of both sets would indicate CP- violation. A non- vanishing value for 
either Q'(6y) or Q'(6y) destroys the hermiticity of the effective theory. 

In the context of SU(2)l ® U(1)y symmetry, the couplings by and by can be realized as first 
order corrections arising from dimension-six operators such as Ff^^F'^"^'^^ or F^^^F'^'^^'^^ where 
$ is the usual Higgs doublet, F^u the field strength tensor and F^i, its dual [51]. Of course, 
higher-order terms may also contribute. Equivalently, the relevant coupling constants may be 



^ Terms not respecting current conservation make vanishing contributions once a gauge boson couples to light 
fermions, as at least one of them must in realistic experimental situations. 
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thought of as momentum dependent form factors. However, for a theory with a cut-off scale A 
large compared to the energy scale at which the scattering experiment is to be performed, the 
form-factor behaviour would be very weak and hence can be neglected for our study. Keeping in 
view the purported higher-order nature of the anomalous couplings, we shall retain only terms 
up to the linear order in all our expressions. 

It is worthwhile to note here that, while our Eq. [1] is the most general expression for the 
VVH vertex, consistent with Lorentz invariance and current conservation, the process under 
consideration, namely e~^e~ —>■ Hff, can, in fact, receive anomalous/no n-SM contributions 
from additional possible operators in an effective theory. Examples include contact interactions 
such as the dimension-6 ffVH operator [55] 



could contribute. The second term can arise from ultraviolet-complete theories, such as a theory 
with a Z' and accommodating a Z'Z'H vertex in the limit of a very heavy Z'. The first one, 
on the other hand, would require a Z'ZH vertex as well. Other constructions, such as theories 
living in higher dimensions, could also lead to such terms in an appropriate approximation [54]. 
While the Qej terms can be neglected in an effective theory approach, the Xjr terms obviously 
have to be included in the most general analysis of the process e+e~ ffH [56-58]. Luckily, 
the contributions of such terms, arising say from a Z' exchange, to the e^e~ ffH amplitude, 
have the same structure as that due to some of the terms in our anomalous vertex (Eq. [1]), as 
long as AjF are flavour universal. For a generic theory — say, with a Z' whose couplings to 
fermions are not flavour universal — the two contributions may be distinguished from each other 
by a comparison of possible differences in different channels. In the present work we desist from 
doing so and thus implicitly assume a flavour universality of the underlying UV-completion 
(say, the Z' couplings). The only remaining dimension-6 operator that is relevant to the given 
process is of the form {iD^e^ {D'^(j)), where i and e are fermionic SU{2) doublet and singlet 
respectively. However, owing to a different chirality structure, it does not interfere with the SM 
amplitude for a massless fermion and hence the corresponding contribution is highly suppressed 



Various terms in the effective VVH vertex have definite properties under the discrete trans- 
formations CP and T, where T stands for the pseudo-time reversal transformation, one which 
reverses particle momenta and spins but does not interchange initial and final states. Table [T] 
shows the behaviour under the transformations, CP and T of various operators in the effective 
Lagrangian, involving different coefficients given in the table. 

3 Kinematics and cuts 

In this analysis, we largely consider the case of ILC operating at a center of mass energy of 
500 GeV and focus on the case of an intermediate mass Higgs boson (2mb < mn < 140 GeV), 
for which H ^ bb is the dominant decay mode with a branching fraction > 0.68 [59]. To be 



A2 



where J-'{3 f /e) denotes a SU{2) multiplet. Even dimension-8 terms such as 



9ei 

A4 



($t$) (e^Me) (/^^/) 
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Table 1: Transformation properties of the various operators (identified by their coefficients) in 
the effective Lagrangian. 

specific, we choose the mass of the Higgs boson to be 120 GeV and the 6-tagging efficiency to 
be 0.7. 

To be detectable, each of the final state particles in the process e~^e~ — > ffH{bb), must 
have a minimum energy and a minimum angular deviation from the beam pipe. Moreover, 
to be recognized as different entities, they need to be well separated. On the other hand, if 
the final state contains neutrinos, then the event must be characterized by a minimum missing 
transverse momentum. Quantitatively, the requirements are 



Ef 


> 


10 GeV 


for each visible outgoing fermion 




< 


175° 


for each visible outgoing fermion 


„miss 
Pt 


> 


15 GeV 


for events with z/'s 


ARjj 


> 


0.7 


for each pair of jets 


ARee 


> 


0.2 


for each pair of charged leptons 


ARij 


> 


0.4 


for jet-lepton isolation . 



(3) 



Here (Ai?)^ = (A0)^ + (Ar/)^, A0 and Arj being the separation between the two entities in 
azimuthal angle and rapidity respectively. 

In addition, cuts may be imposed on the invariant mass of the // system to enhance(suppress) 
the contributions coming from s-channel Z-exchange process, namely 



Rl = 

R2 = 



TTlff — Mz 

rUff — Mz 



< 5Tz => select Z-pole , 

(4) 

> 5 =^ de-select Z-pole , 



where Tz is the width of the Z boson. For the uuH final state, the same goal may be reached 
instead by demanding 

-Rl' = Ett < Eff < E%, 

H- H_ H, 

R2' = Eh < Ejj or Eh > E+, 

where E^ = {s + ml - {mz T 5r^)2)/(2v/i). 

As described earlier, given that the anomalous couplings correspond to operators that 
are notionally suppressed by some large scale, we need retain terms only upto the linear order. 
So, any observable O may be expressed as 

Om}) = OsM + E OiBi. 

The possible sensitivity of these observables to Bj, at a given degree of statistical significance /, 
can be obtained by demanding that \0{{Bi}) - Osm\ < f AO. Here O({0}) = Osm is the SM 
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value of O and AO is the fluctuation in the measurement of O, obtained by adding statistical 
and systematic errors in quadrature. For example, for O being the total cross section (cr) or 
some asymmetry [A), we have 

{Aaf = a/C + e^a\ 

1 - 42 ,2 (6) 

{AAf = '-^ + Ui-Ar. 

a L I 

Here, L is the integrated luminosity of the e^e~ collider and e is the fractional systematic error 
in cross section measurements. Throughout our analysis, we shall take / = 3 and e = 0.01. In 
the approximation of retaining only terms upto linear order in anomalous couplings, the SM 
values for a and A are to be used in Eqs. [61 

It may be noted here that in certain cases, such as when 

1 

0"SM > 



£e2 

the fluctuations are dominated by the systematic error and thus 

{Aaf ^ . 

For example, the cross sections with the i?2'-cut for neutrino final state that are used to 
constrain WWH couplings in Sec. [5] satisfy this condition and hence the bounds on these 
couplings are dominated by systematic errors. 



4 Anomalous ZZH Couplings 

The analysis of Ref. [34] for the case of unpolarised beams had revealed that various asymmetries 
probing the ZZH anomalous couplings were in fact proportional to {I'j—r'j) where If {rj) denote 
the coupling of a Z boson to a left- (right-) handed fermion. Use of polarised beams or detection 
of a r with a specific polarisation in the final state can then avoid the cancellation between 
these two terms and may lead to an enhancement in sensitivity. In this section, we analyse 
various observables that can be constructed with the use of these two quantities. 



4.1 Observables 

Starting from various kinematical quantities Cj constructed as various combinations of the 
different particle momenta and their spins, we define observables as expectation values of 
the signs of Cj, i.e. Oi = (sign(Cj)) (Cj's, i 7^ 1). Each of these observables transform in a 
well-defined manner under C, P and T, and within the aforementioned linear approximation, 
may be used to probe the contribution of a given operator(s) in the effective Lagrangian with 
the same transformation propertied. In fact, the observables listed here are the same ones as 
considered in Ref. [34], but here we use them for a specific polarisation of initial beams and 
final state r's. Needless to say, we concentrate on the cases where use of polarisation affords a 
distinct gain in sensitivity. Table [2] lists some of these observables (cross sections and various 
asymmetries), their transformation properties and the anomalous coupling they may constrain; 
and, below, we give a description of the same: 

^Henceforth, we shall interchangeably use the terminology "transformation properties of anomalous cou- 
plings" for those of the corresponding operator in the effective Lagrangian. 
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Pe -PH 




+ 




+ 




^FB 
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(Pe X Ph) 


■Pf 




+ 
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^UD 


^{k) 


4 


[Pe ■ Ph] * 


[{Pe 


X Ph) ■ Pf] 






+ 






-'^comb 


Q{h) 


5 


[Pe-Pf]* 


[{Pe 


X Ph) ■ Pf] 






(8) 






-'4 comb 


Q{b.)Mbz) 



Table 2: Various possible Ci 's, their transformation properties, the associated observables Oi 
and the anomalous couplings on which they provide information. The symbol ® indicates the 

— * — * 

lack of a definitive transformation property. Here, Pg = Pe- — Pe+, Pf = Pf — Pf o-nd pn is the 
momentum of Higgs boson (to be deduced from the measurement of its decay products). 



1. Oi is nothing but the total cross section as obtained with a specific choice of polarisation 

for the initial beams and/or that for a final state r. As we retain contributions only 
upto the lowest non-trivial order in the anomalous couplings (keeping in view the higher 
dimensional nature of their origin) , the differential cross section can be expressed as 

da ^ [{1 + 2 ^{Aav))daov + '^{^av)daQv + ^{bv)daiv 

v=z,w 

+ ?R:{bv)da2v + '^{bv)dasv + '^{bv)dcr4v] ■ (7) 

where, as in Ref. [34], we have assumed that the Higgs is SM-like and hence 

ay = 1 + Aav (8) 

is close to its SM value. As stated before we choose az to be real, hence Q{Aaz) — and 
we denote ^{Aaz) — Aaz- 

2. O2 is simply the forward-backward(FB) asymmetry with respect to polar angle of the Higgs 

boson, namely 

a{cos9H > 0) - a{cos9H < 0) , , 

^--^'-'^-^ - a(cos^^>0) + .(cos^H<0) 

Since C2 = Pe ■ Ph is odd under CP and even under T transformation, this observable 
would thus be proportional to '^{bz). 

3. O3 is the up-down(UD) asymmetry defined in terms of the momentum of the final state 

fermion / with respect to the iJ-production plane: 

^ a(sin(/) > 0) - o-(sin0 < 0) 

~ (7(sin0 > 0) + (7(sin0 < 0) ■ ^ ' 

As C3 = {Pe X Ph) ■ Pf is odd under both CP and T, one may use this to probe ^(bz). 
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4. C4 = [Pe ■ Ph] * [{Pe X Ph) ■ Pf] is eveii under CP and odd under T and thus expected to be 

sensitive to '^{bz)- The corresponding observable 04 is a particular combination of the 
polar and azimuthal asymmetries (designed to increase sensitivity) and is defined as 

^comb — ; ; ; , U-Lj 

O'FU + O-BD + O-pD + CTBU 

where F, B, U and D refer to the restricted phase space as mentioned above in O2 and 
O3. Thus apu refers to the cross section with Higgs boson restricted to be produced in 
forward hemi-sphere with respect to the initial state electron and the final state fermion 
is produced above the i^-production plane. 

5. is yet another asymmetry derived from a combination of polar and azimuthal distribu- 
tions and is given by: 

_ CTf'U + CTb'D — O'F'D — <^B'U 

^ comb — ; ; ; • l-'-^J 

O'F'U + CTB'D + O'F'D + O'B'U 

Here F' (B') refer to the production of / in forward (backward) hemi-sphere with respect 
to the initial state e~, whereas U and D are the same as defined before. This being both 
P- and T-odd and with no specific C transformation, can be used to probe both Q{bz) 
and ?fi{bz). 

Note that the asymmetries O3, O4 and O5 require charge measurement of the final state particles 
and hence events with light quarks in the final state can not be considered for these observables. 



4.2 Effect of Beam Polarization 

For longitudinally polarized beams, the cross section can be written as 



1 

4 
1 

+ 4 



'l+Ve-){l + Ve+)aRR + {l + Ve-)il-Ve+)cJRL 
;i - n-)(l + rc+)crLR + (1 - n-)(l - V,+ )aLL 



where aipt corresponds to the case of the electron (positron) beams being completely left (right) 
polarized respectively, i.e. , Ve- = —1, Pe+ = +1- ctrr, crl and aiL are defined analogously. 
While the ideal case of complete polarisation is impossible to achieve, values of 80% (60%) 
polarisation for e~(e^) seem possible at the ILC [17]. Taking these to be our default values, we 
denote: 

= a{Ve- = -0.8, Pe+ = 0.6) (13) 

and similarly for other combinations, viz. 0"+'+, (T+'~ and a^'^. We concentrate here on the 
observables discussed in the Sec. 14. II for specific polarization combinations. We would find that 
polarization plays a crucial role in probing the CP-odd ZZH couplings, while the improvement 
in sensitivity in the others is only marginal. 

We quote all our results for an integrated luminosity of 500 fb~^ and a degree of statistical 
significance / = 3 assuming the fractional systematic error to be 1% i.e. e of Eq. [6] to be 0.01. 
Denoting the four possible polarization combinations by 

a : {-,+), b : (+,-), 
c : (-,-), d : (+,+), 



(14) 
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we consider two possible ways of dividing the luminosity amongst these runs, namely 

option (i) : £ = 125 fb~^ for each of (a, b, c, d) 

^ \ 200 fb~^ for (a, 6) (15) 
°P*^°^^")^ ^ = I 50fb-i for(c,d) ■ 

While option (i) is a straightforward choice, option (ii) is better appreciated on realizing that 
polarization combinations (c, d) suppress both SM s-channel processes as also lyW^-fusion. 
Thus, although these combinations maybe useful for certain physics beyond the SM, it is not 
certain whether such modes would find favour for generic search strategies. 

4.2.1 Aaz and 3?(6z) 

It is obvious that the contribution of Aa^ would be identical in form to that within the SM 
and, thus, such a coupling can be probed only through a deviation of the cross sections from 
the SM expectations. As for 3fJ(6z), the fact that this term too conserves each of C, P and T, 
renders all asymmetries insensitive to it. Thus, this coupling too needs to be measured from 
cross sections alone. 

Clearly, just one measurement cannot resolve between these two couplings, and this problem 
was faced by the analysis of Ref. [34] as well. Presumably, with beam polarization being 
available, cross section measurement for a variety of polarization states would offer additional 
information. However, as far as the Bjorken process goes, this dependence is, understandably, 
trivial and identical for both of Aa^ and 3^(6^). This is attested to by the first two rows of 
Table [3], which in fact lists different anomalous contributions to cross sections corresponding to 
different initial state polarisation combinations and different final states. For e^e~ — >■ e^e~H 
though, two diagrams contribute, the usual s-channel one and an additional t-channel one {ZZ 
fusion), with the polarization dependence of the latter being grossly different. To accentuate 
this, we may de-select the Z'-pole (the i?2-cut of Eq. Hj) and the corresponding cross sections 
are displayed in Table [31 For completely polarized beams, LR and RL are the CP-eigen 
states, whereas LL and RR states are CP-conjugate to each other. Hence aiL and a^R receive 
additional contribution from '^{hz) and this is reflected in TableO In addition, this contribution 
is proportional to {Ve- + Ve+) and thus would vanish if the average values of this quantity 
vanishes (as, for example, happens in the unpolarized case). More importantly, the 3?(6z) 
contribution to the ZZ-fusion diagram has an opposite dependence on the product {Ve-Ve+) as 
compared to that of the az piece. This may be exploited to construct an appropriate observable, 
namely 

= [15.1(1 + 2Aa,) + 0.038 ^{bz)] fb ^ ' 

That the contribution of 3(6^) in Eq. [12] vanishes identically is easy to understand. Even 
though the relative weights of the two terms in Eq. [16] can be tuned to reduce the coefficient of 
^{bz) further, it is not really necessary, since both Aaz and 3ft(6z) are expected to arise, say at 
one-loop, and hence would have similar order of magnitude. The large difference in the relative 
weights renders 3fJ(&z) almost irrelevant, making it plausible to constrain Aa^ independent of 
^{bz)- A lack of deviation of Oa from its SM value would give a 3a level limit on Aa^ of the 
form 

r 0.038 for option (i) 

\Aaz\ < S n n^Q f f (■■\ 
0.043 for option (nj. 
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Observable 


Description 


0"OZ 






Ola 


o-"'+(i?l;/i,g) 


23.9 


226 





Olb 


a+'-(i?l;/i,g) 


17.9 


169 





O'la 


(T-'+(i?2;e) 


4.04 


1.46 


0.122 


O'l, 


a+'-(i?2;e) 


2.64 


0.715 


-0.122 


O' 


a-'-(i?2;e) 


3.29 


-1.34 


0.855 




a+'+(/?2;e) 


3.09 


-1.45 


-0.855 



Table 3: Various anomalous contributions (as defined in Eq. [7)j to the cross sections a{Rl; fi, q) 
(for fi^ and light quarks in the final state with Rl-cut) and a{R2;e) (for in the final state 
with R2-cut) for the four polarisation states. The rates are in femtoharns, for ^/s = 500 GeV. 

The smallness of the difference in the two limits is but a consequence of the larger error bars 
resulting from smaller luminosities assigned to polarization combinations {c,d). Note that 
the limits are only marginally different from that obtained with unpolarised beams with C = 
500 fb~^, namely of |Aa^| < 0.04 [34jf|. More importantly, though, unlike the bound of Ref. [34], 
the constraint of Eq. [T7]is independent of 



0.015 


*».. 

s. 






0.01 














0.005 




,\ ^^^"^^ 













0.005 








-0.01 










— — — — — — — — — — — — — — - 




0.015 









-0.04 -0.02 0.02 0.04 

Aaz 



Figure 1: The regions in the Aaz — ^{bz) plane consistent with 3a variations in the observables 
Oa,b,c respectively. The region enclosed by all the three sets reflect the overall constraints. The 
ellipse represents the region corresponding to Ax^ = 11.83 obtained using all polarised cross 
sections listed in Table\^ An integrated luminosity of 125 fb~^ for each of the polarisation state 
i.e. option (i) of Eq. \T^has been used. The limits for option (ii) are very similar. 

A different linear combination of the same observables, namely, 

Ob = Oia + Oi,-Q.Q{0[, + 0[,) 

= [-0.31(1 + 2Aaz) + 4133?(6z)] fb ^ ' 

■^Note that the results quoted here for the unpolarised beams differ from those of the Ref. [34] because of 
the difference in the value of branching fraction oi H ^ bb decay mode used there. We have used the value of 
branching fraction to be 0.68 whereas in Ref. [34] it was taken to be 0.9 
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is equally useful as this enhances the 3?(6z) contribution, while essentially getting rid of Aa^. 
This leads to 

f 0.012 for option (i) , , 

\ 0.018 for option (n) 

virtually independent of Aa^. Finally, using the information from the i?l-cut alone, we have 

= [41.8(1 + 2Aa,) + 395 fb 

leading to a correlated constraint in the Aaz-^{bz) plane as displayed in Fig. [H Of the six 
cross sections of Table [3], one may be used to eliminate 3(6^), leaving behind five constraints 
in this plane. Note that we have already used three linearly independent combinations. We 
may, nonetheless, use all five to construct a x^-test. The resultant Sex ellipse (corresponding 
to Ax^ = 11.83) is also displayed in Fig. [H That this ellipse protrudes slightly beyond the 
set of straight lines is not surprising, for the latter denote the 3a constraint on a particular 
combination of the two variables (with complete disregard for the orthogonal combination), 
while the ellipse gives the corresponding bound on the plane. Furthermore, the size and the 
shape of the ellipse demonstrates that the three combinations identified above represent the 
strongest constraints with very little role played by the remaining two. 

Finally, recollect that, for completely polarized beams, LL and RR states are CP- 
conjugate to each other. Hence the difference of a^L and ajm can be used to probe CP-odd 
couplings. Thus, using 

Od = - = [0.2 (1 + 2Aa,) + 0.11 ^{bz) + 1.71 ^(bz)] fb , (21) 

one obtains 

\^{bz)\ < 0.4 for £ = 125fb-^ (22) 

However, a better constraint can be obtained on this coupling with the use of FB-asymmetry 
with respect to polar angle of Higgs boson and a discussion of this follows. 



4.2.2 ^{bz) and '^{bz) 

Next we focus on the role of beam polarisation in exploring both the real and imaginary parts of 
bz- As discussed earlier, the independent experimental probes for these couphngs, namely the 
asymmetries Ap^ and Aud, are proportional to the quantity (/^ — r^) for the case of unpolarised 
beams [34]. For maximally polarized beams, on the other hand, the squared matrix element 
is either proportional to or to and hence the suppression factor is not so severe even for 
moderate polarization. Since > r^, the cross sections are somewhat larger for the polarization 
combination a = (— , +), and hence the corresponding constraints would turn out to be a little 
stronger. 

Imposing the various kinematical cuts of Eq. [3], along-with the Pl-cut of Eq. H] to select 
Z-pole contribution, the forward-backward (referring to the Higgs polar angle) asymmetry for 
different final states and with different polarisation states of the initial beams, can be written. 
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keeping terms up to linear order in anomalous couplings, as 



Ap'j^{Rl - cut) 



and 



Apj^{Rl - cut) 



0.174 g?(&z) -6.14 '^(bz) 

L48 
-6.07 '^(bz) 

L46 
-92.8 Q(bz) 

22A 

-0.0911 3?(6z) +4.43 %bz) 



1.11 



4.4 Qjbz) 

1.09 
67.2 %bz) 
16.8 



(e+e-H) 

iqqH) 

{e+e-H) 

(qqH) 



(23) 



(24) 



Here, each of the numerical factors denote cross-sections in femtobarns with the denominator 
being the SM cross section. Understandably, the FB asymmetry is identical for the case of the 
Z going into muons or a qq pair, and thus the two channels can be added up to obtain the total 
sensitivity. As for the e^e~H channel, the coupling 3^(6^) makes an appearance on account of 
the interference of the t-channel diagram with the absorptive part of the s-channel SM one. 

The FB-asymmetry with final state /I's and q's get contribution only from '^{bz)- Hence we 
define the following observable 



-W.S^ibz) (25) 



which leads to 



\^^(h )\ < f ^-^^^ option (i) 

'^^^'-[ 0.0089 for option (ii). 



(26) 



In Fig. [2], the vertical lines represent the above bounds for option (i). 

Similarly, the up-down asymmetries, with respect to azimuthal angle of final state fermions, 
are given by 



Ajj'^ {Rl - cut) 



and 



-1.43 g?(&z) -0.286 '^(bz) 

-1.49 ^{bz) 
L46 

' 1.12 3fJ(6z) - 0.161 %bz) 



Apj^ (^1 - cut) 



1.11 



1.08 ^(bz) 
L09 



{e+e-H) 
{e+e~H) 



(27) 



f28) 



Since the determination of this asymmetry requires charge measurement of the final state 
particles, we do not consider it for quarks in the final states. Once again, '^{bz) makes an 
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appearance for the e~^e H case on account of the aforementioned interference. Combining both 
polarisation states (— , +) and (+, — ) for final state muons we construct a observable, namely, 

Oud[RI] /i) = A^^^{Rl- /x) - /x) = - 2.01 ^(bz) , (29) 

we may constrain 

Clearly, one expects a nontrivial effect of the beam poalrisation only for asymmetries con- 
structed with the i?l-cut. For the sake of completeness, we also present the up-down asymme- 
tries with the -R2-cut (de-selecting the Z-pole) for the eeH final state, namely 



A^^(i?2;e) 

Aj}b(i?2;e) 
A+^(i?2;e) 



4.3 af?(6z) +0-227 ^(&z) 
404 

3 3ft(6z) -0.227 $5(6z) 
2^64 ' 

4:.Qm(bz) + 1.59 S(&z) 
3^29 

3.82af?(&z) - 1.59 S(6z' 
3l)9 



(31) 



(32) 



With LL and RR initial states being CP-conjugate to each other, it is understandable that 
AuD receives additional contributioiu from ^{hz)- Defining 

OuD{R2]e) = 2A^^(P2;e) + A+b(^2;e) + Af}b(^2;e) + A+^(P2;e) 

= 5.72 3fJ(6z) - 0.005 3(6^) 

one may get rid of this contribution, and thereby constrain 

|5^(MI<|nn^! (33) 
' ^ - I 0.074 for option (n) ^ ' 

at the 3(7 level. The bounds for option (i) are represented by horizontal lines in Fig. [2l 

It is instructive to compare the above sensitivities (see Table Hj) with those possible with 
unpolarised beams [34]. As one can see, for asymmetries with i?l-cut, the enhancement of 
sensitivity to both the CP-odd couplings 9ft(6z) and ^{pz) on using longitudinally polarized 
beams for option (ii) (option (i)) of Eqs. [I5]) is nearly by a factor of 7(5-6), a feat unachievable 
without polarisation. This improvement is indeed due to the circumvention of the vanishingly 
small vector coupling of electron to the Z boson. Our results are compatible with the more 
stringent limits of Ref . [33] , when we remove the effect of kinematical cuts as well as of the use 
of the hb final state and finite 6-tagging efficiency, implemented in our analysis. Of course, a 
similar enhancement for 9^(6^) may also be possible even with unpolarised beams, when Ajjj:, 
for the eeH final state is considered alongwith the P2-cut. It is, nonetheless, interesting to have 
two different measurements measure the same coupling to the same accuracy. Of course, this 
in addition to the big enhancement in sensitivity to the CP-odd couplings 3?(6z) and ^{hz) as 
mentioned above; this, in fact, is not achievable with the Pl-cut. 



*This is analogous to the appearance of Q{bz) in the total cross sections in Sec l4.2.1l 
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Using Polarized Beams 


Unpolarized States 


Coupling 


Limits 


Observable 
used 


Limits 


Observable 
used 


( iFtf ion (11 


( )T~tf 10T1 (111 


1 CYA / 7 1 ^ 

|3fJ(6z| < 


0.067 


0.074 




0.067 




\Si{hz)\ < 


0.17 


0.13 




0.91 






0.011 


0.0089 




0.064 





Table 4: Limits on anomalous ZZH couplings from various observables at 3a level for both 
polarized and unpolarized beams. While an integrated luminosity of 500 fb~^ is assumed for the 
unpolarized run, for the polarized case the same is divided according to the options of Eq. [73 

0.08 
0.06 
0.04 
0.02 



^ -0.02 
-0.04 
-0.06 
-0.08 

-0.015 -0.01 -0.005 0.005 0.01 0.015 

Im(bz) 




Figure 2: The regions in the Q{bz)^{bz) plane consistent with 3cr variations in the observables 
Ofb{R^', yW, q) and Oud{R'2'] e). The vertical and horizontal lines are for respective variation in 
them. The ellipse represents the region corresponding to a 3cr bound in the plane (Ax^ = 11.83^ 
obtained using all asymmetries listed in Eqs. 
division has been assumed. 



'2l\ \2B, EH Option (i) of luminosity 



A further improvement is possible if all the asymmetries listed in Eqs. ( !23p4l [271128]) are 
considered alongwith the three independent combinations obtained from Eqs. [SUby eliminating 
^ipz) therefrom. Once again, a can be constructed leading to a 3cr contour as displayed in 
Fig. [21 

And, finally, we use the remaining information available from Eqs. [31] to probe 3(6^). 
Defining 

0'uD{R'^]e) = A-jj^o{R2-e)- A+^(i?2; e) = - 0.017 3?(6z) + 0.998 , (34) 
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we may impose 

I'^ih)] < 0.22 for £ = 125 fb~^ (35) 

This limit is only marginally different than the one obtained with unpolarised beams using the 
combined asymmetry Acomb of Sec. 14.11 with i?l-cut. We shall see in the next section that a 
better constraint may be obtained on this coupling from the combined asymmetry if the final 
state r helicity is measured. 

4.3 Use of Final State r Polarization with Unpolarised Beams. 

In this section we report on the use of selecting final states with r's in a given helicity state. A 
similar idea was employed in the optimal variable analysis of Ref. [31]. A detailed measurement 
of the decay pion energy distribution [60] or a simpler measurement of the inclusive single 
pion spectrum [61], can yield information on r polarisation. Discussions exist in literature 
on how this can be utilized to sharpen search strategies for charged Higgs boson [60-66] or 
super symmetric partners [67] and even for the measurement of SUSY parameters [68,69]. In a 
similar vein, one can also construct observables using the final state r polarisation to probe ZZH 
couplings. We construct asymmetries for r with definite polarisation which can be measured in 
simple counting experiments and which can catch the essence of the above mentioned optimal 
observable analysis. 

As discussed in Ref. [34], for unpolarized initial and final states, the combined polar- 
azimuthal asymmetry Acomb is proportional to {il + r'^){r'j — £j) and the up-down asymmetry, 
A(7d(0) is proportional to {i^ — r^^rj — ij-). Thus, for leptonic final states, both these asym- 
metries suffer a suppression (this is particularly important for both Aud and Acomb which are 
impossible to measure with hadronic final states). Hence, the measurement of the final state r 
polarisation would lead to an enhancement in these symmetries with a consequent improvement 
in the sensitivity to both of the T-odd anomalous couplings, namely 9(&z) and 5J(6z). Further, 
since i"^ > r^, one gets a slightly higher gain in sensitivity with final state r in a negative 
hehcity state. 

To demonstrate this, we construct various asymmetries (listed in Table [2]) for both left- and 
right-handed r in the final state. These are the same as described in Sec. 14.11 but defined for a 
particular helicity of r rather than for a specific initial beam polarisation state (for this analysis 
we take the initial state to be unpolarised). Again, after imposing the kinematical cuts of Eq. [3] 
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alongwith the -Rl-cut, these read 



A' 



-0.527 ^{hz) 
0.495 ' 
0.388 ^{hz) 
0.365 ' 
-127 %hz) 

0.495 
1-01 %bz) 
0.365 ' 
1.18 53(6^) - 



(36) 



0.3 ^{hz) 



0.495 

-0.868 %hz) - 0.221 ^{hz) 



0.365 



where the numbers represent the partial (total) cross sections in femtobarns. As expected, the 
asymmetries are, generically, in the opposite sense for left- and right-handed r's (and, hence, 
would have largely cancelled each other). Thus, being able to isolate r's of a given hehcity 
indeed can add to the sensitivity and the slightly larger value of means that the final state 
with negative helicity r's will be even better. Rather than effect a detailed analysis, we make a 
simplifying and illustrative assumption that final state r~'s in a specific (negative or positive) 
helicity state can be isolated with very high purity at the cost of an efficiency of 40% or even 
20%. Under this assumption, we list, in Table [5|, the best limits that are possible for the T-odd 
couplings ^{pz) and 3?(6^) using asymmetries of Eq. [36] for negative helicity r's. (The positive 
helicity r's would lead to very similar (marginally weaker) bounds, and we omit these for ease 
of presentation.) Corresponding limits obtainable with the r final state but without using 
knowledge of r helicity are also shown for comparison. 





Using Pol. of final state r 


Unpolarised r's 


Coupling 


Limits 


Observable 


Limits 


Observable 


40% eff. 


20% eff. 




0.11 


0.15 


comb 


0.35 


A comb 




0.28 


0.40 




0.91 


AuD 



Table 5: Limits on anomalous ZZH couplings from various observables at 3a level with nega- 
tively polarised and unpolarised t's (i.e. without using knowledge of t helicity), with Rl-cut. 



It is clear from Table [5] that the limits on Q{bz) and ^{bz) can be improved by the mea- 
surement of the final state r polarization. Preliminary results of our analysis were presented 
in Ref. [70]. A similar observation had been made earlier in the context of optimal observable 
analysis [31]. It is to be noted here that the unpolarized measurements with ecH final state, 
with i?2-cut gives a better sensitivity to 3?(6z) as discussed earlier. It is, however, nice to have 
more than one observable determining a given coupling. We would like to stress here that even 
with an efficiency of isolating events with a given r~ helicity as low as 20%, this method affords 
an increase in the sensitivity for '^{bz) by as much as a factor ~ 2. 
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Figure 3: 3a blind regions in the 3^(6^) — Q{bz) for a r-helicity isolation efficiency of 40%. 
and an integrated luminosity of 500 fh^^ . Panel s. (h) is for left- (right- )handed t~ 's. The 
horizontal, vertical and oblique lines correspond to Acomh, ^ud aiT'd A'comb respectively. The 
ellipse combines all three and represents the blind region in the plane. The inner ellipse of 
panel (c) combines both sets, and the outer represents the constraints without using r-helicity 
information. 

Fig. [3] displays the region in 3ft(6z) — ^{bz) plane that can be probed using the asymmetries 
of Eq. [36] for r~'s in a specific helicity state, for an isolation efficiency of 40%. Once again, 
a test can be constructed trivially. Note that the constraint from A'comb are in opposite 
sense for left- and right-handed r's, resulting in a rotation of the two 3a ellipses with respect to 
each other. One may then combine the information obtained from the two final states and the 
results thereof are shown in the last panel. The fact of the two individual ellipses being only 
slightly rotated with respect to the 3?(6z) axes means that this exercise of combining the two 
leads to only a moderate improvement. It is interesting, though, to compare the result with the 
corresponding region when the r-helicity information is unavailable. The difference is obvious. 

4.4 Use of final state r polarisation for polarised beams 

Having established that each of beam polarization and measurement of r-helicity can lead to 
substantial improvements, the natural question relates to the combination of the two effects. 
Recall that the up-down asymmetry Aud is proportional to {l^ — r'^) (/^ — r^) while the combined 
asymmetry Acomb is proportional to (/g -|- r^) (fl — r^). Thus, while the isolation of events with 
a final state r in a definite helicity state would enhance both, the use of polarised beams will 
only enhance Aud- Hence in this section, we concentrate on the improvement in sensitivity to 
'^{pz) (which is probed by Aud)- 

The up-down(UD) asymmetry for r's in negative helicity state, with negatively polarized 
e~ and positively polarized e"*" beams (i.e. the polarisation state 'a' as mentioned in Sec. 14.20 
is given by 

A-^i(Rl-r,) = Z^l^iMi). (37) 

UDK , LJ Qggg 

Once again, assuming that r^^'s can be isolated with an efficiency of 40% (20%), the associated 3 
a limits on ^{bz) from this single measurement alone for an integrated luminosity L = 200 fb~^ 



18 



(i.e. using only part of the data available for option (ii) of the luminosity divide) reads 



~ [ 0.054 for 40% efficiency 

- 1 0.077 for 20% efficiency. ^^^^ 

Comparing this to Table HI the large gain in sensitivity is obvious. And that too with only part 
of the data. This is significantly better than the maximal sensitivity available for unpolarized 
beams, namely \^{bz) \ < 0.067 obtainable for the e~^e~H final state with i?2-cut and isolating r- 
helicities with a 40% efficiency. Once the data for the other combinations of beam polarizations 
and the r-helicity, viz. 



4.1 ^{bz) 

0.627 
4.17 ^{bz) 



A-ui{Rl;rn) = (39) 



-3.02 lk{bz) 



0.462 

are taken into account (using a test), one obtains, for a 40% isolation efficiency, 

,^ ~ f 0.032 for option (i) , , 

[ 0.040 for option (n) , 

with the numbers for an isolation efficiency of 20% being a little worse. The combined use of 
beam polarisation and r-helicity measurement thus plays a very productive role. 

To summarize the entire section, we have demonstrated that the CP-odd and T-odd ZZH 
couplings can be probed far better by the use of polarised beams and/or the information of final 
state polarisation. The sensitivity to CP-even and T-even anomalous ZZH couplings, (Aa^ 
and '^{bz))i on the other hand, show only a marginal improvement in their sensitivity limits. 
However, use of polarised beams helps to constrain these couplings independent of each other. 



5 Beam Polarization and the WWH Couplings 

A study of the anomalous WWH couplings is possible via the process e^e~ vvH. However, 
since it receives contributions from both the t-channel WW fusion diagram and the s-channel 
Bjorken diagram, this determination needs knowledge of the anomalous ZZH couplings, which 
fortunately, can be measured well from measurements of other final states. It may be argued 
that, for completely polarized beams, the cross section aiR gets contribution from both the 
Bjorken and fusion diagrams, whereas only the first contributes to (Trl, and hence it should be 
possible to use cross sections with different polarisation combinations to reduce the contam- 
ination due to the anomalous ZZH couplings. It should be noted however, that in realistic 
situations one would not have 100% beam polarisation, and thus this effect cannot be entirely 
neutralized. 

Apart from this possible contamination, the determination of WWH anomalous coupling 
suffers from one more limitation. The presence of a pair of neutrinos deprives us of the full 
knowledge of the momenta of the final state fermions and, thus, does not allow construction of 
T-odd observables. Total rates and forward-backward asymmetry with respect to polar angle of 
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the Higgs boson (each for different combinations of beam polarisations) are the only observables 
available in the present case. Using the notation of Eq. [T4]for the polarisation states, the states 
'c' and 'ci' do not contribute to the process e"*" e~ — > wH, for the case of 100% polarisation. 
Even with the assumed values of 80% and 60% polarisation for the e~ and beams respectively, 
these two combinations will correspond to rather small cross sections. Therefore, we restrict 
ourselves to the data accrued from the other two polarisation combinations viz. a and h. On 
imposition of the RV- or the i?2'-cut, the cross sections are now given bj|^ 

(T-'+(i?l'; z/) = [9.09 + 17.6 Aa^ + 0.60 3?(Aaiy) + 83.8 - 3.62 

- 0.48 9(AaM/) + 0.90 ^{hw) + 1-51 "^{hw)] fb 

(T+'-(i?l» = [6.63 + 13.2 Aaz + 0.02 3?(AaH/)+ 63.4 3?(6z)- 0.10 $5(6z) 

- 0.01 $5(AaH/) + 0.03 ^{hw) + 0.04 '^{hw)] fb (41) 

a-^+{R2';u) = [102 - 0.45 Aa^ - 17.6 K(6z) - 0.31 $5(6^) 
+ 205 ^{Aaw) - 38.3 ^(bw)] fb 

(T+'-(i?2';z/) = [3.23 + 0.78 Aaz + 4.31 3ft(6z) + 5.68 sR(AavF) - 1.06 3fJ(6vi/)]fb 

That the WWH couplings make a small contribution for the (+, — ) case is easy to under- 
stand. Furthermore, the contributions from '^{bw,z) are small as these have to be proportional 
to the absorptive part of the Z-propagator in Bjorken diagram. As can be expected, the sensi- 
tivity to the WWH couplings is enhanced if one can successfully eliminate the contribution of 
the Z-diagram (this has the further advantage of eliminating the z/^ and z/^ events). However, 
since even the -R2'-cut (see Eq. |5]) cannot eliminate the Z diagram entirely, as is seen from 
Eq. UH we consider a particular linear combination of the cross sections, namely 

O2UA = a"'+(i?2» +4 a+'-(i?2» 

= [r]i + 115 + 228 ^{Aaw) - 42.5 ^{bw)] fb, (42) 

r/i = 2.66 Aaz- 0.36 3?(6z) -0.31 $5(6z). 

Thus the contamination from ZZH couplings is contained entirely in rji. Using the analysis 
of the previous section (which did not involve the uuH final state), we have |?7i| < 0.13(0.14) 
respectively for options (i) and (ii) of the luminosity division (Eqs. [T5i) . In other words, the 
uncertainty due to a lack of precise knowledge of the ZZH couplings is reduced to negligible 
proportions enabling us to constrain a combination of ^{Aaw) and 3?(&i^) virtually independent 
of ZZH couplings: 

|2 ^{Aaw) - 0.37 ^(bw)] < 0.040 (0.036) (43) 

for the two choices of luminosity division among different polarisation modes concerned (see 
Fig. H]). The relatively small difference between the sensitivities indicates that the luminosity 
division is not very crucial as long as a sufficiently large fraction is devolved into the canonical 
choices of (+, — ) and (— , +). It may be also noted here that due to large value of SM cross 
section, the fluctuation in the cross section with _R2'-cut given by Eq. [6] is dominated by the 
systematic error as mentioned in Sec. [3l Thus small changes in luminosity or cross section is 



^ As stated earlier in Sec. O we may choose az to be real, without any loss of generality. But once we make 
this choice, all the WWH couplings (including Aaiy) may be complex. 
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not going to make much difference to the sensitivity hmits of the couphngs that are probed by 
cross section a~'^{R2'] u) or 02va- 

While it is true that we can only probe a combination of these two couplings using 02va^ 
use of beam polarization clearly helps to reduce contamination from ZZH couplings to this 
determination. It is to be noted that since only a particular linear combination of 9ft(Aaiy) and 
'^{pw) can be probed, limits on them are strongly correlated. Of course, it is possible to derive 
a constraint on the orthogonal combination from the data of Eq. but it is too weak to be of 
any relevance. If we make a further assumption of only one of these couplings being non-zero, 
we would obtain 

|3fJ(AavF)| < 0.020(0.018) for 3?(6h/) = , . 

\^{hw)\ < 0.110(0.095) for 3fJ(Aaiy) = ^ ' 

with the two limits corresponding to options (i) and (ii) of Eqs. [El Although a comparison with 
the results obtained earlier [34] shows only a marginal improvement in the individual limits, the 
determination is now free of uncertainty coming from contamination from the ZZH couplings. 
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Figure 4: The region in the ^{Aaw) — ^{bw) plane consistent with 3a variations in 02ua /^^^ 
an integrated luminosity of 125 fb~^. The vertical line shows the sensitivity limit on Aaz from 
Ftg.m 

At this stage, it is intriguing to consider the consequences of an additional assumption 
(made in Ref. [34]) of Aa^y = Aa^, which is found to be true in some cases [71] and would 
be motivated by SU{2) x f/(l) invariance of the effective theory. By itself, this would impose 
a bound on Aaw, courtesy Eq. [T7] and hence that would lead to a closed shown in 

Fig. m However, one should also note that SU{2) x f/(l) invariance would further equate hw 
and hz, and that the correlation between constraints in the 3?(6vi/)^3^(Aaty) plane (Fig. H]) is 
in the opposite sense to that in the 3?(6^)-Aa^ plane (Fig. [1]). Thus, an assumption of such an 
invariance would lead to far stronger constraints. In particular, the improvement is dramatic 
for 3ft(6vK). 

As for the other CP-even WWH couplings, namely ^^(Aavi/) and ^{pw), note that their 
contributions arise from the interference of the lyW^-fusion diagram with the absorptive part 
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of the Z-propagator in the Bjorken diagram. Hence the corresponding terms appear only in 
total rate with i?l'-cut and being proportional to the width of Z-boson, are very small. Since 
the presence of two neutrinos in the final state does not allow us to construct any T-odd 
observables, this study is virtually insensitive to these two couplings. Of course, an assumption 
of SU{2) X U{1) invariance would change matters drastically. 

What remains is to investigate bw, and this being a CP-odd coupling, can be probed through 
the forward-backward (FB) asymmetry with respect to the polar angle of Higgs boson. These 
asymmetries, RV- or and i?2'-cuts can be expressed as 



^FB 




^) 


1 

9.09 


-2.29 ^(bz) - 36.9 '^(bz) 


+ 0.57 3fJ(6vK) - 0.47 '^(bw)] , 


-^FB 


[RV- 


^) 


1 

6.63 


-0.064 + 27.1 %bz) 


+ 0.02 3?(&H/) -0.01 55(6vf) 


4-'+ 
^FB 


(i?2'; 


^) 


1 r 

102 ^ 


5.17 $5(6z) + 7.83 ^{bw) 






^FB 


(i?2'; 


^) 


1 

3.23 


2.1 53(6^) + 0.22 53(6vk) 







(45) 



where the numbers once again represent the corresponding cross sections in femtobarns. With 
the last two of Eqs. HS] involving just a single WWH coupling, these can be combined to 
eliminate the remaining dependence on the ZZH anomalous couplings to leave us 

A,^^^ ^ 13 A]^^+(i?2»-A+^-(i?2» 
= 0.009 + 0.93 5>(6vf)- 

The contribution to A2mix from '^{bz) (which, incidentally, can be probed very accurately, see 
Table |4]) may be neglected, leading to 

|5>(&vk) I < 0.50(0.44) (47) 

for C = 125 (200) fb^^ for each of the (+, — ) and (— ,+) polarization combinations. While 
it may seem that the improvement is marginal when compared to the sensitivity achievable 
with unpolarised beams — |53(6vk)| ^ 0.46 for a total luminosity of 500 fb~^ [34] — note that, 
unlike the older analysis, the current sensitivity is independent of any other coupling. In other 
words, the use of beam polarisation has allowed construction of an observable that can isolate 
the contribution of $5(6^/). 

The anomalous coupling ^{bw) being a T-odd coupling has no contribution to FB asymme- 
tries with i?2'-cut and only a small contribution to FB asymmetries with i?l'-cut through the 
interference of the lyW^-fusion diagram with the absorptive part of Bjorken diagram. Hence 
these asymmetries are not a good probe of 3^(6^/). Thus the process under consideration can 
not be used to probe any of the T-odd couplings in the WWH vertex. 

To summarize the results of this section: use of beam polarisation allows us to obtain 
limits on T-even couplings ^{Aaw), ^{bw) and '^{bw) independent of the ZZH couplings, 
these being listed in Table El Using linear combinations of our observables corresponding to 
different polarisation combinations for initial states, the contamination from ZZH coupling can 
be reduced to negligible amount, something that was not possible with the unpolarised beams. 
Further, even though this measurement uses only two of the combinations, an equal division of 
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Using Polarized Beams 


Unpolarized States 


Coupling 


Limits (for given C) 


Observables 
used 


Limits 


Observables 
used 


zUU ID 


izo ID 


|3?(A«n-)| < 


0.018 


0.020 




0.019 


a""i"'\B2': p) 


mbw)\ < 


0.095 


0.11 




0.10 






0.44 


0.50 


A2 ■ 


0.40 





Table 6: 3 cr limits on anomalous WWH couplings from various observables with polarized and 
unpolarized beams. 



the total luminosity among all four already gives optimal results. The limits on 3?(Aavi/) and 
"Slipw) are highly correlated whereas '^{pw) is constrained independent of any other coupling. 
There is no observable to constrain the T-odd couplings and thus we conclude that this process 
is not a good probe for these couphngs. Hence one has to look to other processes to probe these 
couplings. For example, this problem may be overcome in the 67 collider as shown in Rcf. [72]. 
The WWH couplings are not contaminated by the ZZH couplings in the process studied there 
viz. (e~7 — > vW' H). Using this process in conjunction with the PT conjugate process, the 
authors were able to construct observables that depend only on one of the anomalous WWH 
couphngs and hence the hmits obtained on each of the WWH couphngs were independent of 
the others. The anomalous VVH couplings at 67 collider have been also studied in Rcf. [73]. 
Finally, an assumption of SU{2) x f/(l) invariance would drastically improve the constraints 
for virtually all the couplings, whether they be WWH or ZZH. 



6 Dependence on ^/s and inclusion of ISR/Beamstrahlung 
Effects 

All our analysis so far has been performed for a fixed centre of mass energy, namely \/i = 
500 GeV. Clearly, the total cross section is a function of energy, and the functional form 
may depend on the presence of (and the identity of) any anomalous coupling, owing to their 
higher- dimensional nature. Thus, the sensitivities could, in principle, depend on the choice 
of \fs. Furthermore, for the processes involving both s- and t-channel diagrams, the relative 
importance of these two parts of the amplitude is, in fact, energy dependent, and the generic 
enhancement of the f-channel cross section with increasing beam energies may, in principle, lead 
to an improvement in the sensitivity to those couplings primarily constrained using observables 
with i?2-(i?2')-cut. 

Even for a nominally fixed beam energy, the ^/s available to an individual hard scatter- 
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ing event is generically less than this value owing to the ubiquitous initial state radiation 
(ISR) — which is nothing but the bremsstrahlung radiation by the incoming particles — or beam- 
strahlung, which is a name for the radiation from the beam particles due to its interaction with 
the (strong) electromagnetic fields caused by the dense bunches of the opposite charge in a col- 
lider environment. Consequently, it is important to investigate the possibly detrimental effects 
of such eventualities. 




Figure 5: ^/s-variation of (aoz(-R2; e) (the SM part in the cross section a{R2] e) = cr(e+e^ 
e'^e~H) with R2-cut), (Tovy(-R2'; z/) and aiw{R2'] v) (the coefficients of '^{/^aw) and 3fJ(6iy) 
respectively in the cross section for e^e~ —>■ vvH with R2'-cut). 



To this end, we begin by studying the y^-dependence of the observables used in this paper. 
For simplicity, we shall restrict ourselves, in this section, to unpolarised scattering, with the 
results for polarised beams expected to be similar. Moreover, we shall concentrate on observ- 
ables defined with the i?2-(i?2')-cut. For example, the ZZH coupling Aaz is best probed by 
the total cross section with i?2-cut for electron final state i.e. a{R2; e) = a{e~^e'' e^e~H) 
with i?2-cut. Fig. \5\ shows the variation of the SM part of this cross section (ctoz of Eq. [7]) 
with v^|§ Recall that Aaz simply rescales the SM part of the cross section (Eq. [7]) and hence 
the sensitivity is determined simply by the corresponding number of events, namely it scales as 

1 /2 

Nqz' . With aoz{R2; e) having a maximum at y/s ^1.1 TeV, this would be the optimal beam 
energy to probe Aaz, with the maximal improvement (compared to the results quoted above) 
being by about 40%. 

Fig. also shows the y^-variation of the SM and ^{bw) contributions (i.e. aow and aiw 
respectively) to the cross section a{R2'; u) = a{e^e^ vvH) with /22'-cut which is the best 
probe for the WWH couplings "^{Aaw) and 3?(6vk). While the monotonic increase in aiw would 
seem to suggest that increasing ^/s would readily lead to an improvement in the sensitivity to 
"^{hw), note that this increase saturates and, furthermore, that aow increases at least as fast. 
Consequently, for moderate changes in y^, any improvement or otherwise is expected to be 
marginal at best. 

^Note that ctoz has been scaled by a factor of 50 to fit in the same figure. 
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Moving to asymmetries, the up-down asymmetry Au£){R2; e) of the final state fermion 
with respect to the if-production plane, and the forward-backward asymmetry with respect 
to polar angle of the Higgs boson, Afb{,R'^'] ^) have been used to constrain ^{bz) and '^{bw) 
respectively [34]. Expressing these as 

AuDiR2;e) = A«^3fJ(6z) 

AFBiR2';iy) = 53(6^) + S(M (48) 
the coefficients A^^, and Ap^ are plotted in Fig. [6] as a function of ^/s. The figure clearly shows 
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Figure 6: y/s-variation of the coefficients A^^ in up-down asymmetry and the Ap^ in forward- 
backward asymmetry as defined in Eq. 48 of the text. Both plots are with R2-(R2')-cut (de- 
selecting Z pole) imposed. 

that the sensitivity to 3^(6^) is expected to improve at higher y/s, while, ^/s = 500 GeV is an 
optimal choice for '^{bw) and higher energies would only tend to deteriorate the sensitivity. 
The arguments above are refiected by Table [7] which summarizes the sensitivity limits at the 
3 a level, with an integrated luminosity of 500 fb~^, for different center of mass energies. 

The above analysis did not take into account either of ISR and beamstrahlung. We now 
proceed to do so, using the structure function formalism [48] to incorporate these effects. The 
differential scattering cross section for a given process e~ (pi) + {p2) — * -^(7) can be expressed 
as: 

c?(T[e+e" X(7)] = fe/e{xi) fe/e{x2) c?a[e+e" X]{s) , 

where the electron luminosity function fe/e{x) describes the probability of finding an electron 
with a momentum fraction x of the nominal beam energy, or, in other words, the probability 
with which an electron energy E = A/i/2 emits one or more photons with total energy (1 — x)E 
resulting in the reduction of its energy to E^. = xE. Hence the square of the effective cm. 
energy s can be expressed as: s ~ X1X2S. 

Using the Weiszacker- Williams approximation, one can write down the luminosity function 
for ISR as [47] 



Je/e ' 



X 



16 



(8 + 3/5)(l 



4(1 + x) 



(49) 



25 



Coupling 




3 a limit at 

^/s = 0.5 

TeV 


3 a limit at 

Vs = 1 
TeV 


3 (7 limit at 

Vs = 3 
TeV 


Observable 
used 


|Aaz| 


< 


0.040 
0.043 


0.030 
0.031 


0.049 
0.039 


a{R2; e) 
^ISR+Beam.(^2;e) 


mbz)\ 


< 


0.067 
0.075 


0.028 
0.032 


0.015 
0.018 


AuD{R'2;e) 
^ISR+Beam.^^2;e) 


mAaw)\ 


< 


0.019 
0.019 


0.016 
0.017 


0.016 
0.016 


a{R2'; v) 
^ISR+Beam.^^2» 


mbw)\ 


< 


0.10 
0.11 


0.082 
0.084 


0.084 
0.083 


^ISR+Beam.(^2» 


l^(MI 


< 


0.40 
0.43 


0.42 
0.41 


0.89 
0.71 


ApB{R2'\y) 
^BR+Beam.(^2» 



Table 7: Individual sensitivity limits (assuming only the relevant coupling to he non-zero) at the 
3(7 level for an integrated luminosity of 500 fb^^ on various anomalous couplings at different 
cm. energies without and with the inclusion of ISR and beamstrahlung effects. The results 
correspond to unpolarised beams. 



where 

/^^^flog^-l) , (50) 

and «em is the fine-structure constant. 

The beamstrahlung spectrum depends on electron beam energy E, and the parameters such 
as the number of electrons per bunch iVg, the bunch dimensions (for a Gaussian bunch profile) in 
both the longitudinal direction (cr^) as well as in the transverse directions {(Ty,x)- K is convenient 
to introduce a "beamstrahlung parameter" T [49] given by: 

where r^, is the classical electron radius, and rrif. its mass. The electron spectrum f^^/'^"^{x) 
describing the effects of beamstrahlung, can be written in a closed analytical form for T < 10 
[49]. Armed with all these, the expression for the electron spectrum function, including both 
ISR and beamstrahlung effects, may be written as [48] : 
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where flf^iO is as given in Eq. HQ] above, whereas /^g^™ ^ x) is that given by Eq. 22 of 
Ref. [49]. 

We now analyse how the ISR and beamstrahlung effects modify the sensitivity hmits of 
^{bz), for which we had observed an improvement in sensitivity at higher ^/s as hsted in 
Table [3 In our analysis we use the beamstrahlung parameters to be [74] 

(Tz = 30 /im for all energies, 

T = 0.3, 1.0 and 8.1 for E^^ = 0.5, 1.0, and 3.0 TeV respectively. (52) 
Fig. [7] compares the total rates with and without ISR effect for different final states. The 
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Figure 7: y/s-variation of cross sections with and without the ISR effects for final states with 
muons or electrons. The figure also shows the points at -y/i = 0.5 TeV and 1.0 TeV when 
beamstrahlung effects are also included. Beamstrahlung effects are negligible at y/s = 0.5 TeV. 

effective cm. energy s of the electron beam decreases due to ISR. As a result, away from 
the threshold, the ISR effects increase the rates of pure s-channel processes, such as e"*" e~ 
fi~^ fi~ H . Near the threshold, however, the cross section rises with the effective ^/s and hence 
the ISR effects cause a decrease in the rate. For the process under discussion, the cross-over 
from one behaviour to the other takes place at about 300 GeV. On the other hand, {e^e~H)- 
production receives both s- and t- channel contributions, and hence there is no such cross-over 
with ISR effects always decreasing the rates. We thus expect these effects to reduce somewhat 
the improvement observed above in probing ^{bz) at higher energies. 

At a/s = 500 GeV, beamstrahlung effects are negligible whereas ISR effects cause 10% to 
15% change in cross section (see Fig. [7]). However, this does not cause significant changes in 
sensitivities as both the SM and anomalous parts are affected in much the same way. At higher 
cm. energies both ISR and beamstrahlung effects are nontrivial and significant. However, once 
again, the effect is in the same direction, and roughly of the same magnitude, for the SM and in 
the presence of anomalous couplings (Fig. [7j). Although the figure demonstrates this for total 
cross sections, the story is similar even for partial cross sections and for the other couplings as 
well. 
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Since the constraints on ^{bz) are the only ones to improve significantly with increasing ^/s 
(see Table [7]), we concentrate on the effects of ISR and beamstrahlung on this coupling. With 
^{bz) being best probed by the up-down asymmetry (with i?2-cut) for electrons in {e~^e~H) 
production, the ISR and beamstrahlung effects can be summarised as listed in Table [HI 





UD-asymmetry(y4[/o(i?2; e)) 


No ISR & Beam. 


With ISR 


With ISR & Beam. 


0.5 TeV 


1.16 ^{bz) 


1.13 ^{bz) 


1.1 


1 TeV 


2.00 ^{bz) 


1.94 ^(bz) 


1.85 ^{bz) 


3 TeV 


6.29 ^(bz) 


5.60 ^(bz) 


4.23 ^(bz) 



Table 8: Up-down asymmetry for final state electron for R2-cut(Amj{R2; e) ) with and without 
ISR and beamstrahlung effects at different ^/s's. 

As is readily seen, the effects are almost negligible even for y/s = 1 TeV, and only marginally 
important for ^/s = 3 TeV. The consequent shift in the sensitivities are summarised in Table 
[71 The smallness of the effects can be understood by realizing that, on the imposition of the 
i?2-cut (de-selecting the Z-pole), the energy dependence of the cross sections (total or partial) 
is only logarithmic. Further, both SM and anomalous parts have a similar dependence. Hence 
although the cross sections are affected significantly there is little effect on the sensitivity limits 
of the anomalous parts after inclusion of ISR and beamstrahlung effects. 

We conclude this section by making a few general observations : 

• With increasing energy, the observables with the i?2-cut imposed gain more in sensitivity 
as compared to those defined with the -Rl-cut. 

• Observables with i?l-cut (selecting Z-pole) are affected more by ISR and beamstrahlung 
corrections because of the usual s-channel suppression, whereas the observables with R2- 
cut (de-selecting Z-pole) have only logarithmic A/s-dependence and hence do not suffer 
as significant corrections. 

• By using higher cm. energies, we gain significantly in the sensitivity to 3^(6^) (upto a 
factor of 2 at y/s = 1 TeV as compared to the case of 500 GeV). Neither the ISR nor the 
beamstrahlung effects change the sensitivity significantly. 

• There is no significant gain in the sensitivity to anomalous WWH couplings even at very 
high energies. 

• In totality, we do not gain much in the sensitivity by going to higher cm. energies. Thus, 
running the collider at lower energies (say y/s = 500 GeV), but with polarised beams is 
more beneficial for the study of anomalous VVH couplings. 
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7 Summary 



With the plethora of the gauge boson couphngs with a Higgs that an effective theory allows, 
it is always going to be a comphcated task to unravel these. It was shown in Ref. [34] that 
it is possible to construct observables with definite CP and T transformation properties that 
are sensitive to a single anomalous coupling (the one with matching CP and T properties) and 
hence these observables may be used to obtain robust constraints on the couplings in a model 
independent way. With polarised beams, we have the possibility of constructing even more 
such observables. 

We have performed our analysis by imposing kinematical cuts on different final state particles 
to reduce backgrounds. In addition, we also include the reduction (in the event rates) caused 
by considering only the events wherein the H decays into a hh pair (with branching fraction 
~ 0.68), with a realistic 6-tagging efficiency of 70%. This renders, to a great extent, our 
estimates of sensitivity quite realistic. 

Almost all of the asymmetries pertaining to the probe of the ZZH vertex, that were con- 
structed in Ref. [34] for unpolarised beams and neglecting the polarisation of the final state 
particles, are proportional to the difference between the squared right and left handed coupling 
of the Z to charged leptons. Consequently, one expects the sensitivity to the corresponding 
anomalous couplings to be enhanced once one uses a specific polarisation, thereby overcoming 
this cancellation between the (almost) equal left and right handed couphngs of the Z boson to 
leptons. Indeed, we observe this. For example, the use of beam polarization alone leads on the 
one hand, to the disentangling of the couplings Aa^ and 3?(6^) (not possible for unpolarised 
beams), and on the other, to an improvement in the sensitivities to CP-odd ZZH couplings 
by a factor of upto 5-7. 

With the helicity of a r in the final state decipherable from the momentum distribution of 
the charged prongs in its decay, one can construct asymmetries pertaining to r's with a specific 
helicity. Assuming that it would be possible to isolate events with r's in a definite helicity 
state with an efficiency of (say) 40%, we use such asymmetries to demonstrate that sensitivity 
to '^{hz) and "^{hz) could be improved by a factor of about 3. While an earlier optimal 
observable analysis [31] had indeed investigated the use of r helicity, our analysis differs in that 
we have constructed an observable that can be measured in a simple counting experiment and 
catches the essence of the optimal observable. Also worth noting is that the use of final state r 
measurement along with polarised beams, allows us to improve on the sensitivity for 3?(6z) by 
a factor of about 2. 

As far as the constraints on the anomalous WWH couphngs are concerned, the problem of 
a smaller number of observables due to the presence of missing neutrinos in the the i'DH{W) 
still exists. However, by constructing new observables that are combinations of the observables 
corresponding to different polarisation states, the contamination from the anomalous ZZH 
couplings is highly reduced in case of T-even WWH couplings which is a great advantage 
over the use of unpolarised beams. For T-odd WWH couplings though, we conclude that this 
process is not a good probe for them because of non- availability of any T-odd observable. Hence 
one has to look to other processes to probe these couplings. 

We have also investigated possible gains in sensitivity by going to higher centre of mass 
energies. This though renders ISR and the beamstrahlung effects progressively more important. 
As a matter of fact, even for = 500 GeV, the ISR effects can affect cross sections with the 
i?l-cut by 10-15%. Fortunately, the ISR and beamstrahlung effects on sensitivity limits is very 
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modest, as both the standard model and the anomalous contributions change by similar amounts 
causing only small numerical effects, especially on the asymmetries. As for the variation in 
sensitivity with ^/s ranging from 300 GeV to 3 TeV, the observables with the i?l-cut (selecting 
the Z-pole) are naturally associated with a larger cross section (and, hence, better limits) 
for lower energies. Thus, a machine operating at ^/s = 350 GeV would do even better for 
certain couplings than the present case. This is in agreement with the findings of Ref. [34]. 
For observables constructed with the i?2-cut, however, one can obtain an improvement in the 
sensitivity limit (say) for ^{bz) by upto a factor of 2 at cm. energy 1 TeV compared to the 
ones possible at 500 GeV [34]. 

In conclusion, we find that use of beam polarisation and final state r polarisation can result 
in significant advances in the search for anomalous ZZH couplings at the ILC. For WWH 
couplings, though there is no significant change in the sensitivity limits, the T-even ones can 
be now probed virtually independent of the anomalous ZZH couplings. The inclusion of ISR 
and beamstrahlung effects, changes the individual cross sections but has very little effect on 
the sensitivities. Further, going to higher beam energies, leads to only modest improvements. 
Thus, as far as this sector of the theory is concerned, there is a strong case for use of beam 
polarisation and measurement of final state fermion polarisations wherever possible, but there 
is no real gain in going to higher energies. 
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Note: During the last stages of finalising our manuscript, a paper dealing with related issues 
containing a very detailed analysis [75] on probing the coefficients of the CP-even, dimension 6 
operators in the anomalous VVH vertices has appeared. 

References 

[1] P.W. Higgs, Phys. Rev. Lett. 13 (1964) 508; Phys. Rev. 145 (1966) 1156; F. Englert and 
R. Brout, Phys. Rev. Lett. 13 (1964) 321; G.S. Guralnik, C.R. Hagen and T. Kibble, Phys. 
Rev. Lett. 13 (1964) 585. 

[2] J. Gunion, H. Haber, G. Kane, and S. Dawson, The Higgs Hunter's Guide, Addison-Wesley, 
Reading (USA), 1990. 

[3] M. Gomez-Bock et al, J. Phys. Conf. Ser. 18, 74 (2005) | arXiv:hep-ph/0509077| . 



30 



[4] A. Djouadi, Phys. Kept . 457, 1 (2008) |arXiv:hep-ph/ 0503172] ; Phys. Kept. 459, 1 (2008) 
I arXiv:hep-ph/0503173| . 



[5] For a review, see for example, R. M. Godbole, [arXiv : hep-ph /0205114 ], Volume 4, Jubilee 
Issue of the Indian Journal of Physics, pp. 44-83, 2004, Guest Editors: A. Raychaudhuri 
and P. Mitra. 

[6] ATLAS Collaboration, Technical Design Report, CERN-LHCC-99-14 and CERN-LHCC- 
99-15; CMS Collaboration, Technical Design Report, CMS-LHCC-2006-21. 

[7] ALEPH, DELPHI, L3 and OPAL Collaborations, The LEP Working 
Group for Higgs Boson Searches, Phys. Lett. B565, 61 (2003); see also 



|http://lephiggs.web.cern.ch/LEPHIGGS/ 



[8] LEP Electroweak Working Group, 'http://lepewwg.web.cern.ch/LEPEWWG/ J. Alcaraz 
et al. [LEP Collaborations], CERN-PH-EP/2007-039, arX iv:0712.0929l [hep-exj. 

[9] R. Akers et al. [OPAL Collaboration], Z. Phys. C61, 19 (1994); D. Buskulic et al. [ALEPH 
Collaboration], Z. Phys. C62, 539 (1994); M. Acciarri et al. [L3 Collaboration], Z. Phys. 
C62, 551 (1994); P. Abreu et al. [DELPHI Collaboration], Z. Phys. C67, 69 (1995); 
D. Antreasyan et al. [Crystal Ball Collaboration], Phys. Lett. B251, 204 (1990); P. Franzini 
et al, Phys. Rev. D35, 2883 (1987); J. Kalinowski and M. Krawczyk, Phys. Lett. B361, 66 
(1995) |arXiv:hep-ph/9506291|; D. Choudhury and M. Krawczyk, Phys. Rev. D55, 2774 
(1997) [arXiv:hep-ph/9607271] . 

[10] M. Carena, J. R. Ellis, A. Pilaftsis and C. E. M. Wagner, Nucl. Phys. B586, 92 (2000). 

[11] A. Mendez and A. Pomarol, Phys. Lett. B279, 98 (1992); J. F. Gunion, H. E. Haber 
and J. Wudka, Phys. Rev. D43, 904 (1991); J. F. Gunion, B. Grzadkowski, H. E. Haber 
and J. Kalinowski, Phys. Rev. Lett. 79, 982 (1997) [hep-ph/9704410|; I. F. Ginzburg, 
M. Krawczyk and P. Osland, ,hep-ph/0211371j I. F. Ginzburg and M. Krawczyk, Phys. 
Rev. D 72, 115013 (2005). 

[12] See e.g. M. Drees, R.M. Godbole and P. Roy, Theory and phenomenology of sparticles, 
World Scientific, 2005; H. Baer and X. Tata, "Weak scale Super symmetry: From superfields 
to scattering events," Cambridge, UK: Univ. Pr. (2006) . 

[13] G. Abbiendi et al. [OPAL Collaboration], Eur. Phys. J. C 37, 49 (2004) 
|arXi v:hep-ex/04060 57|; ALEPH, DELPHI, L3 and OPAL Collaborations, The LEP 
Working Group for Higgs Boson Searches, LHWG-Note 2004-01, CERN-PH-EP-2006- 
001; S. Schael et al. [ALEPH Collaboration], Eur. Phys. J. C 47 (2006) 547 



|arXiv:hep-ex/0602042| . 
[14] E. Accomando et al... Report of 'Workshop on CP studies and non-standard Higgs physics. 



CERN-2006-009, arXiv:hep-ph/0608079, 



[15] See for example, J. R. Forshaw, Pramana 63, 1119 (2004) and discussion therein; see 
also D. Ch oudhury, T. M. P. Tait and C. E. M. Wagner, Phys. Rev. D65, 053002 (2002) 
|arXiv:hep-ph/0109097|. 



31 



[16] E. Accomando et ai, Phys. Kept. 299 (1998) 1 
Saavedra et ai, |hep-ph / 1 63 1 5 ; T. Abe et ai, 



|hep-ph/0109166 



arXiv:hep-ph/9705442|; J. Aguilar- 



hep-ex/0 106055 -58; K. Abe et al, 



[ir 

[18 
[19 

[2o; 

[21 

[22 
[23 
[24 

[25; 



[26; 
[27; 
[28; 

[29; 

[30 
[31 
[32; 



J. Brau et al, SLAC-R-857; A. Djouadi et a/.. larXiv:0709.1893l [hep-ph]. 

G. Weiglein et al. [LHC/LC Study Group], Phys. Kept. 426, 47 (2006) 
I arXiv:hep-ph/0410364| . 



For a review, see, e.g., R. M. Godbole et al, arXiv:hep-ph/0404024 and references therein. 

V. D. Barger et ai, Phys. Rev. D 49 (1994) 79 |arXiv:h ep-ph/9306270|, A. Skjold and 
P. Osland, Phys. Lett. B 311 (1993) 261 |a rXiv:hep-ph/9303294j ; Nud. Phys. B 453 
(1995) 3 [arXiv:hep-ph/9502283] . 

A. Djouadi and B. Kniehl, Correlations in Higgs production and decay as a probe of CP- 
violation in the scalar sector, pubhshed in the Proceedings of the Workshop on "e"'"e~ 
Colhsions" , Report DESY-93-123C. 

D. J. Miller et al, Phys. Lett. B 505 (2001) 149 |arXiv:hep-ph/0 102023]. 
M. T. Dova, P. Garcia-Abia and W. Lohmann, |arXiv:hep-ph/0302113 



P. S. Bhupal Dev et ai, Phys. Rev. Lett. 100, 051801 (2008) [a rXiv:0707.2878l [hep-ph]]: 
[arXiv:0710.2669 [hep-ph]. 

R. M. Godbole and P. Roy, Phys. Rev. Lett. 50 (1983) 717; P. Kalyniak, J. N. Ng, P. Za- 
karauskas, Phys. Rev. D 29, 502 (1984), J. C. Romao and A. Barroso, Phys. Lett. B 
185, 195 (1987). W. Ki lian, M. Kramer and P. M. Zerwas, Phys. Lett. B 373, 135 (1996) 
jarXiv : hep-ph /9512355| and referneces therein. 

K. Hagiwara and M. L. Stong, Z. Phys. C 62 (1994) 99 [a rXiv:hep-ph/9309248| . 
R. Rattazzi, Z. Phys. C 40 (1988) 605. 

H. J. He, Y. P. Kuang, C. P. Yuan and B. Zhang, Phys. Lett. B 554, 64 (2003) 
|arXiv:hep-ph/02Tl 229| . 

K. Hagiwara, S. Ishihara, R. Szalapski and D. Zeppenfeld, Phys. Rev. D 48, 2182 
(1993); K. Hagiwara, R. Szalapski and D. Zeppenfeld, Phys. Lett. B 318, 155 (1993) 
|arXiv:hep-ph/9308347| . 

G . J. Gounaris, F. M . Renard and N. D. Vlachos, Nucl. Phys. B 459, 51 (1996) 
|arXiv:hep-ph/9509316|. 

K. Hagiwara, S. Ishihara, J. Kamoshita and B. A. Kniehl, Eur. Phys. J. C 14, 457 (2000) 
I a rXiv:hep-ph/0002043 | . 

D. Chang, W. Y. Keung and I. Phillips, Phys. Rev. D 48 (1993) 3225 
[arXiv:hep-ph/9303226j . 



32 



T. Han and J. Jiang, Phys. Rev. D 63, 096007 (2001) [arXiv:hep- ph/001127 1|. 

S. S. Biswal, R. M. Godbole, R. K. Singh and D. Choudhury, Phys. Rev. D 73, 035001 
(2006) [Erratum-ibid. D 74, 039904 (2006)] |arXiv:hep-ph/0509070|. 

S. S. Biswal, D. Choudhury, R. M. Godbole and R. K. Singh, Pramana 69, 777 (2007). 

D. Atwood and A. Soni, Phys. Rev. D 45, 2405 (1992). 

S. Y. Choi, D. J. Miller, M . M. Muhlleitner and P. M. Zerwas, Phys. Lett. B 553 (2003) 
61 |arXiv:hep-ph/02 10077). 

C. P. Buszello, I. Fleck, P. Marquard and J. J. van der Bij, Eur. Phys. J. C 32 (2004) 
209 |arXiv:hep-ph/0212396|; C. P. Buszello and P. Marquard, in Ref. [14]; C. P. Buszello, 



P. Marquard and J. J. van der Bij, |arXiv:hep-ph/0406181 
M. Bluj, Acta Phys. Polon. B 38, 739 (2007). 

R. M. Godbole, D. J. Miller, S. Moretti and M. M. Miihlleitner, Pramana 67 (2006) 617; 
zbzd. in Ref. [14]. 

B. Zhang, Y. P. Kuan g, H. J. He and C. P. Yuan, Phys. Rev. D 67, 114024 (2003) 
|arXiv:hep-ph/0303048| . 

S. M. Lietti, A. A. Natale, C. G. Roldao and R. Rosenfeld, Phys. Lett. B 497 (2001) 243 
|arXi v:hep-p h/0009289| . 

R. M. Godbole, D. J. Miller and M. M. Muhlleitner, JHEP 0712, 031 (2007) 
[a rXiv:0708.0458l [hep-ph]]. 

T. Plehn, D. L. Rainwater and D. Zeppenfeld, Phys. Rev. Lett. 88, 051801 (2002) 
I a rXiv:hep-ph/0105325 ] . 

C. P. Buszello and P. Marquard, |arXiv:hep-ph/0603209[ 

P. Niezurawski, A. F. Z arnecki and M. Krawczyk, Acta Phys. Polon. B 36, 833 (2005) 
|arXiv:hep-ph/0410291| . 

E. A. Kuraev and V. S. Fadin, Sov. J. Nucl. Phys. 41, 466 (1985) [Yad. Fiz. 41, 733 
(1985)]. 



M. Drees and R. M. Godbole, Z. Phys. C 59, 591 (1993) [arXiv:hep-ph/9203219] . 
P. Chen, Phys. Rev. D 46, 1186 (1992). 

B. A. Kniehl, Nucl. Phys. B 352, 1 (1991); Nucl. Phys. B 357, 439 (1991) 

M. C. Gonzalez-Garcia, Int. J. Mod. Phys. A 1 4, 3121 (1999) |arXiv :hep-ph/9902321j ; 
V. Barger et al, Phys. Rev. D 67, 115001 (2003) )arXiv:hep-ph/0301097[^ 

[52] J. F. Gunion, H. E. Haber and J. Wudka, Phys. Rev. D 43 (1991) 904. 



33 



[53] A. Mendez and A. Pomarol, Phys. Lett. B 272, 313 (1991). 

[54] D. Choudhury, A. Datta and K. Huitu, Nucl. Phys. B 673, 385 (2003) 
I a rXiv:hep-ph/0302141 ] . 

[55] W. Buchmuller and D. Wyler, Nucl. Phys. B 268, 621 (1986). 

[56] B. Grzadkowski and J. Wudka, Phys. Lett. B 364, 49 (1995) |arXiv:hep-ph/9502415| . 

[57] W. Kihan, M. K ramer and P. M. Zerwas, Phys. Lett. B 381, 243 (1996) 
I arXiv:hep-ph/9603409| . 

[58] K. Rao and S. D. Rindani, Phys. Lett. B 642, 85 (2006) |arXiv:hep-ph/0605298| ; Phys. 
Rev. D 77, 015009 (2008) [arXiv:0709.2591 [hep-ph]] 

[59] A. Djouadi, J. Kahnowski and M. Spira, Comput. Phys. Commun. 108, 56 (1998) 
I arXiv:hep-ph/9704448] . 

[60] K. Hagiwara, A. D. Martin and D. Zeppenfeld, Phys. Lett. B 235, 198 (1990). 

[61] D. P. Roy, Phys. Lett. B 277 (1992) 183. 

[62] B. K. Bullock, K. Hagiwara and A. D. Martin, Phys. Rev. Lett. 67 (1991) 3055. 

[63] B. K. Bullock, K. Hagiwara and A. D. Martin, Nucl. Phys. B 395 (1993) 499. 

[64] S. Raychaudhuri and D. P. Roy, Phys. Rev. D 52 (1995) 1556. 

[65] S. Raychaudhuri and D. P. Roy, Phys. Rev. D 53 (1996) 4902. 

[66] G. Bella, Nucl. Phys. Proc. Suppl. 40 (1995) 475. 

[67] M. Guchait and D. P. Roy, arXiv:0808.0438l [hep-ph]. 

[68] M. M. Nojiri, Phys. Rev. D51 (1995) 6281. 

[69] R. M. Godbol e, M. Guchait and D. P. Roy, Phys. Lett. B 618, 193 (2005) 
[arXi v:hep-ph/0411306| . 

[70] S. S. Biswal, D. Choudhury, R. M. Godbole and Mamta. [a?Xiv:0710.2735l [hep-ph]. 

[71] A. Pilaftsis and C. E. M. Wagner, Nucl. Phys. B553, 3 (1999); M. Carena, H. E. Haber, 
H. E. Logan and S. Mrenna, Phys. Rev. D 65, 055005 (2002). 

[72] D. Choudhury and Mamta, Phys. Rev. D74 (2006) 115019 | arXiv:hep-ph/0608293l ; Pra- 
mana 69, 795 (2007). 

[73] I. Sahin, Phys. Rev. D 77, 115010 (2008) [ arXiv:0802.0293l [hep-ph]]. 

[74] R. W. Assmann et ai, SLAC-REPRINT-2000-096. 

[75] S. Dutta, K. Hagiwara, and Yu Matsumoto. arXiv:0808.0477y l [hep-ph]. 



34 



